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Abstract. Action at distance in Newtonian physics is replaced by 
finite propagation speeds in classical post-Newtonian physics. As 
a result, the differential equations of motion in Newtonian physics 
are replaced by functional differential equations, where the delay 
associated with the hnite propagation speed is taken into account. 
Newtonian equations of motion, with post-Newtonian corrections, 
are often used to approximate the functional differential equations. 
In |H] a simple atomic model based on a functional differential 
equation which reproduces the quantized Bohr atomic model was 
presented. The unique assumption was that the electrodynamic 
interaction has finite propagation speed. Are the finite propagation 
speeds also the origin of the gravitational quantization? In this 
work a simple gravitational model based on a functional differential 
equation gives an explanation of the modified Titius-Bode law. 

1. Introduction 

In the last two centuries several attempts were made to express and 
explain the distribution of the planetary orbits and other relevant quan- 
tities using integer numbers. Titius (1772) and Bode (1776) (see for 
instance O 1221) proposed the law describing the mean distances of 
planets from the Sun of the general form 

r n = a + bc n , 

where r n means distance characterized by an integer number n. The 
constants a, b and c have no convincing physical meaning, neither have 
the empirical correlations with definite parameters for a given system. 
Therefore, this law has raised many discussions. Nevertheless, it played 
a positive role not only in predicting unknown planets, but also in 
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stimulating many researches to further work in this direction. In |TT] 
and [2Uj there are some reviews about the attempts to test and to 
explain the Titius-Bode law. 

As a precursor of the modified Titius-Bode law, Gulak |H] proposed 
that the orbital distances are given by r n = (n + l/2)r or r n = nr , 
where r$ is a characteristic of a given system. Here n needs not to in- 
crease by 1 in going from one planet or satellite to another one. Gulak 
[TUj found a theoretical support to his previous results by constructing 
an equation of the Schrodinger type. In this way, he tried to introduce 
the macroquantization of orbits in a gravitational field. 

In the last years the idea that of a quantization of the gravitational 
field has been constated. In words of Halton Arp (see [2]): An unex- 
pected property of astronomical objects (and therefore an ignored and 
suppressed subject) is that their properties are quantized, for instance 
the redshifts of galaxies. The most astonishing result was then pointed 
to by Jess Artem, that the same quantization ratio that appeared in 
quasar redshifts appeared in the orbital parameters of the planets in 
the solar system. Shortly, afterward Oliveira Neto et al. [16], Agnese 
and Festa 0, L. Nottale et al. [HJ [T5] and A. and J. Rubcic |2U1 I2T] 
independently in Brazil, Italy, France and Croatia began pointing out 
similarities to the Bohr atom in the orbital placement of the planets. 
Different variations of the Bohr-like r n = n 2 or r n = n 2 + n/2 fit the 
planetary semimajor axes extremely well with rather low "quantum" 
numbers n. It is clear that the properties of the planets are not ran- 
dom and that they are in some way connected to quantum mechanical 
parameters both of which are connected to cosmological properties. 

Action at distance in Newtonian physics is replaced by finite prop- 
agation speeds in classical post-Newtonian physics. As a result, the 
differential equations of motion in Newtonian physics are replaced by 
functional differential equations, where the delay associated with the 
finite propagation speed is taken into account. Newtonian equations 
of motion, with post-Newtonian corrections, are often used to approx- 
imate the functional differential equations, see, for instance, [13 IH El 
El HI UH1 EH]- In [HI a simple atomic model based on a functional dif- 
ferential equation which reproduces the quantized Bohr atomic model 
was presented. The unique assumption was that the electrodynamic 
interaction has finite propagation speed, which is a consequence of the 
Relativity theory. An straightforward consequence of the theory devel- 
oped in jH], and taking into account that gravitational interaction has 
also a finite propagation speed, is that the same model is applicable to 
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the gravitational 2-body problem. In the following section we present a 
simple gravitational model based on a functional differential equation 
which gives an explanation of the modified Titius-Bode law. 

2. The retarded gravitational 2-body problem 

We consider two particles of masses m and M, with m <C M, in- 
teracting through the retarded inverse square force. The force on the 
mass m exerted by the mass M is given by 




Figure 1 . The retarded gravitational 2-body problem. 

The force acts in the direction of the 3-vector r, along which the mass 
M is "last seen" by the mass m. The 3-vector r may be represented 
by 

r = r M (t - t) - r m (t), 

where i\m(£) and r m (t) denote respectively the instantaneous position 
vectors of the mass M and the mass to, respectively, at time t, and r 
is the delay, so that — t) is the "last seen" position of the mass 
M. Assuming that the two bodies are in rigid rotation with constant 
angular velocity u, and referring back to Fig. 1, we have, in 3-vector 
notation, 

r m = r i [cos ujt i + sin Ult j] , 

and 

v M = — r 2 [cosuj(t — r) 1 + smu(t — r) j]. 
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Hence, the 3-vector r is given by 

r = [— r 2 cos u(t — r) — r\ cos uot] T + [— r 2 sin u(t — r) — r\ sin uot] j, 

Now, we introduce the polar coordinates (r, 9) and define the unitary 
vectors 1 = cos 9 i+sin 9 j and n = — sin 9 i+cos 9 j. By straightforward 
calculations it is easy to see that the components of the force (Q) in the 
polar coordinates are 

Mm Mm 
r r = G — — r ■ 1 = (— r 2 cosfwr) - ri)G — — 

and 

Mm . Mm 

re = G — 5— r • n = r 2 sin (cur) G — 5— 

The equations of the movement are 

(2) mr — mr9 = F r , 

(3) mr9 + 2mr9 = Fg. 

The second equation (jHJ) can be written in the form 

(4) -— = - — (mr 2 ^) = F e = r 2 sm^G^— 
r at r at r d 

If we accurately study equation JH) we see that the analytic function 
sin(cijr) has a numerable number of zeros given by 

(5) LOT = klT , 

with k G Z, which are stationary orbits of the system of equations (J2J) 
and (J3J). When cur 7^ A;7r we have a torque which conduces the mass 
m to the stationary orbits without torque, that is, with ujt = kir. In 
fact the stationary orbits are limit cycles in the sense of the qualitative 
theory developed by Poincare, see [T7] . 

This is a new form of treating the gravitational 2-body problem from 
a dynamic point of view instead of from a static point of view, as it 
has been made up to now. Moreover, in this model the delay r is not 
small, in fact 

kir kn kT 
T ~ ~ ~ 2vr/T ~ ~Y ' 

where T is the time taken by the mass m to complete its orbit, i.e., 
the period of revolution. Therefore, the delay is a multiple of the half- 
period T/2. 

On the other hand, in a first approximation, the delay r can be equal 
to r/c (the time that the field uses to goes from the mass M to the 
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mass m at the speed of the light). In this case, from equation (JIjJ) we 
have 

6 r = — = - . 

uj c 

Taking into account that u = v$/r, from © we have Vg/c = kn. 
However, from the Relativity theory we know that Vg/c < 1, then we 
must introduce a new constant g in the delay. Hence, r = gr/c and 
the new equation (jUJ) is 

kn or 

7 r = — = U — , 

UJ c 

and now Vg/c = kn/g, i.e. vg = knc/g and from (JJJ) we also have 
r = kirc/(guj). In our model case of a classical rigid rotation we have 
9 = ujt with uj > 0. Therefore, 9 = uj and 9 = 0. Hence, equation (j3J) 
for cg>t = A;7r is 

2mfuj = 0, 

which implies f = and r = where is a constant for each k. On 
the other hand, equation (J2J) for cur = kit takes the form: 

8 -mrfl 2 = -m-^ = -r 2 (-1 n - r^G—r- « -r G — 5—, 

assuming that r ~ ri due to r 2 <C ri in the case that m <C M. 

From the definition of angular momentum L = mr 2 9 = mr 2 uj = 
mrvg we have that vg = L/(mr). Substituting this value of Vg into 
equation (jSJ) we obtain r = L 2 /(GMm 2 ). The energy of the mass m 
(substituting the values of vg and r) is given by 

mvj Mm G 2 M 2 m 3 

( 9 ) E = ^ G = T72 — • 

2 r 2L Z 

The angular momentum for ujt = kn is 

kirc L 2 

(10) L = mvgr = m- 



g GM m 2 ' 

which is an equation for the angular momentum. Isolating the value 
of L we obtain L = (GM m g) / (kirc) . If we introduce this value of the 
angular momentum in the expression of the energy Q we have 

(11) E~=gk. 
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The adimensional constant g cannot be unambiguously determined, 
and it must be computed according with the experimental data. How- 
ever, one may speculate by considering the similarity between the grav- 
itational constant force F g and the electrodynamic force F e . The ab- 
solute ratio of the forces is 

F g _ GMm 



F e e 2 /47T£o' 

where Eq is the electric permittivity constant of vacuum. By introduc- 
ing the well-known fine structure constant a defined by 

e 2 

a 



4ne hc 

the ratio may be expressed by 

F g _ 1 (GMm\ _ a g 
F e a \ he ) a ' 

where a g is the adimensional gravitational fine structure constant, see 
[T5j . However, this value of the adimensional gravitational fine struc- 
ture constant does not agree with the experimental data given in the 
works [TJ HH EH I20l ED • If we recall the expression of the energy 
levels for the electrostatic interaction given by Bohr in 1913 

E = --m(ac) — , 
2 n 1 

a straightforward generalization gives us that the expression of the 
energy levels for the gravitational interaction 

(12) E = -~m(a g c) 2 ± 

Therefore, comparing with ()11|) (identifying n = \k\) we must impose 
that the constant g = k 2 gi where g\ = n/a g and then the energy takes 
the form (|12|). Therefore we have found the value of the adimensional 
constant g\ and consequently the expression of the delay r which is 

or , 9 qir , 9 iir 

(13) T = — = k 2 — = k 2 . 

c c a g c 

In fact, we generalize the expression of the delay for the electrodynamic 
interaction given in [S]. From the found value of the angular momentum 
and the value of vg = a g c/ k we have 

GMmg GMmk a„c 

(14) L = — = = mvgr = m— — r. 

kne a g c k 
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Isolating the value of r from equation (fHJ) and identifying n = \k\, we 
arrive to the radii of the stationary orbits 



(15) 



GMn 2 

9 9 ' 

dig & 



which depends on M. Equation (JT5j) agrees with the experimental 
data given in the works El EH 1201 E] where a 9 is empirically 
determined. One of the most important differences between the electro- 
dynamic interaction and the gravitational interaction is that the unit 
of mass starts from zero but the basic unit of charge never changes. For 
this reason it was difficult to detect the gravitational quantization. The 
dependence of the stationary orbits on M and the possible dependence 
of a g on the masses M and m makes it difficult to find a unique law for 
the distribution of planets in the solar system. This is the reason that 
implies the using, in some models, different expressions for the interior 
planets and the exterior planets, see [2H1 EE] If all the planets had 
the same mass m it would be easy to find a general law of distribution 
similar to the quantization given in atomic models. 

As a consequence of the values of vg and r n we have that the period 
T of revolution of the planets is proportional to n 3 because 

_ 2vrr _ 27r ^p^ _ 2nGMk 3 2tcG M n 3 

J- n 7Z 7- 



v e ^ a 3 g c 3 o^c 3 

Therefore, as it must happen, it is satisfied the third Kepler's law, i.e., 
the ratio r 3 /T% does not depend on n. 

Summarizing, with the found delay definition (|13|). the model pre- 
sented in this work explains the modified Titius-Bode law faithfully. 
The gravitational quantization is, in fact, the first approximation in 
the value v/c of the delay in the gravitational interaction. This first 
approach to the macroquantization of orbits is confirmed by the ob- 
served data analyzed in Oliveira Neto et al. [IB] . Agnese and Festa 
[I], L. Nottale et al. O US] and A. and J. Rubcic [HHH]- In these 
works different values of a g are found. For instance, in (201 121] the 
value of the gravitational fine structure constant is a g = a/(2irf), 
where a is the electro dynamic fine structure constant, and / is an adi- 
mensional constant determined for each planet. In Q the value of the 
gravitational fine structure constant determined according with the ex- 
perimental data is a g = 1/2086, which explains the distribution of all 
the planets in the solar system. An open important problem is the 
determination of a g in function of both masses M and m and other 
universal constants. 
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3. Concluding remarks 

In |Sj the atomic Bohr model is completely described by means of 
functional differential equations. It is important to stand out that what 
we will carry out in the following section is not a post-Newtonian ap- 
proach in which r is small. This is what has been made up to now and 
in the mentioned works [3J HJ E] . In this work, we will accept that 
the laws governing the movement have a delay (a delay that does not 
need to be small) and we will find a solution of the functional differen- 
tial equation in a very simple case. In this work we have obtained the 
Newtonian approximation of the gravitational field taking into account 
that the gravitational interaction has finite propagation speed, which is 
a consequence of the Relativity theory. This Newtonian approximation 
is a reminiscent of quantization of the gravitational field. The quanti- 
zation of the gravitational field must be obtained using the Einstein's 
field equation and the delay, which must appear in a natural way in 
this equation. 
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